/•' 


► 


WL-TR-9  3-4060 


AD-A265  906 


ANALYTICAL  MODELS  FOR  DESCRIBING  TH 
MECHANICAL  BEHAVIOR  OF  FIBER  REINFORCED 


DR.  JAMES  M.  WHITNEY 


UNIVERSAL  TECHNOLOGY  CORP . 
4031  COLONEL  GLENN  HWY 
DAYTON  OH  45431-1600 


SEPTEMBER  1992 

FINAL  REPORT  FOR  06/01/92-09/30/92 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  IS  UNLIMITED. 


MATERIALS  DIRECTORATE 

WRIGHT  LABORATORY 

AIR  FORCE  MATERIEL  COMMAND 

WRIGHT  PATTERSON  AFB  OH  45433-7734 


93-13638 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  in  connection  with  a  definitely  Government-related 
procurement,  the  United  States  Government  incurs  no  responsibility  or  any 
obligation  whatsoever.  The  fact  that  the  government  may  have  formulated  or 
in  any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not 
to  be  regarded  by  implication,  or  otherwise  in  any  manner  construed,  as 
licensing  the  holder,  or  any  other  person  or  corporation;  or  as  conveying 
any  rights  or  permission  to  manufacture,  use,  or  sell  any  patented  invention 
that  may  in  any  way  be  related  thereto. 


This  report  is  releasable  to  the  National  Technical  Information  Service 
(NTIS).  At  NTIS,  it  will  be  available  to  the  general  public,  including 
foreign  nations. 


This  technical  report  has  been  reviewed  and  is  approved  for  publica 

tion. 


1  >  ft  oSUc- 

RICHARD  B.  HALL 
Materials  Research  Engineer 
Mechanics  &  Surface  Interactions  Branch 


L.  SCOTT  THEIBERT,  Chief 
Mechanics  &  Surface  Interactions  Branch 
NonmetaJlic  Materials  Division 


ILES  E.  BROWNING,  Chief 
Nonmetallic  Materials  Division 
Materials  Directorate 


If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing 
list,  or  if  the  addressee  is  no  longer  employed  by  your  organization  please 
notify  WL/MLBM  ,  WPAFB,  OH  45433-  7734  to  help  us  maintain  a  current 
mailing  list. 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by 
security  considerations,  contractual  obligations,  or  notice  on  a  specific 
document. 


Form  Approved 
OM8  Nc  0704 -01 88 


REPORT  DOCUMENTATION  PAGE 


Public  reporting  burden  for  this  collection  ot  information  is  estimated  to  average  1  hour  per  response,  including  the  time  tor  reviewing  w 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  tho  collection  of  information  Send  comments  regarding  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Sconces.  Directorate  (or  Informa 
Davis  Highway,  Suite  1204,  Arlington.  VA  22202-4302,  and  to  the  Office  ot  Management  and  Budget.  Paperwork  Reduction  Project  (( 

situations.  searching  existing  data  sources. 

Out  ben  estimate  c?  any  other  aspect  of  rh»s 
bon  Operations  ana  Reports  Jefferson 

)7G4  0188),  Washington.  DC  20603 

1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

April  1993  Final  Jun  92  -  Sep  92 

4.  TITLE  AND  SUBTITLE 

Analytical  Models  for  Describing  the  Mechanical  Behavior  of  Fiber  Reinforced 
Composites 

5.  FUNDING  NUMBERS 

C:  F336 I 5-K9-C-5609 

PE  62102F 

FR  2418 

TA  04 

WU  68 

6.  AUTHOR(S) 

Dr.  James  M.  Whitney 

7.  PERFORMING  ORGANIZATION  NAMES(S)  AND  ADDRESS(ES) 

Universal  Technology  Corp.,  4031  Colonel  Glenn  Hwy.,  Dayton  OH  45431-1600 
Dr.  James  M.  Whitney,  University  of  Dayton,  Dayton  OH  45469-0001 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

WL/MLBM 

2941  P  St  Ste  I 

Wright-Patterson  AFB  OH  45433-7750 

10.  SPONSORING  /  MONITORING 

AGENCY  REPORT  NUMBER 

WL-TR- 93-4060 

11.  SUPPLEMENTARY  NOTES 

Work  performed  by  J.  M.  Whitney  under  subcontract  to  UTC. 

12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (Maximum  200  words) 

This  report  covers  three  topics  in  mechanics  of  advanced  composites:  effect  of  fiber  cross-section  on  stiffness  and  strength 
of  laminated  composites;  failure  criteria  for  composites  under  combined  loading;  and  stress  analysis  models  for  thick 
composites. 

In  the  area  of  fiber  cross-sectional  effect,  micromechanics  analysis  addresses  the  effect  of  elliptical  fibers  on  the  stiffness 
and  in-plane  stress  concentrations.  Numerical  results  show  that  fiber  cross-sectional  geometry  has  only  a  modest  influence 
on  laminate  elastic  constants. 

The  second  topic  addresses  the  interaction  between  longitudinal  compression  and  transverse  loading  (tension  and 
compression).  Data  obtained  with  a  mini-sandwich  beam  are  used  to  generate  failure  stresses.  Classical  lamination  theory-, 
including  thermal  residual  stresses  due  to  cure,  is  utilized  for  determining  ply  stresses.  Experimental  results  reveal  little 
interaction  between  longitudinal  compression  and  transverse  stresses  as  long  as  the  transverse  load  is  below  uniaxial 
transverse  failure. 

Under  the  topic  of  thick  composites,  the  effect  of  shear  deformation  and  transverse  normal  stress  on  the  cylindrical  bending 
of  laminated,  anisotropic  plates  subjected  to  a  uniform  lateral  load  is  investigated.  Closed  form  solutions  indicate  that  the 
effect  of  transverse  normal  stress  on  maximum  deflection  is  small,  even  for  relatively  thick  plates. 

14.  SUBJECT  TERMS 

Composites,  miconechani'-s,  elliptical  fibers,  failure  criteria,  thick  composites,  combined 
loading,  compression,  shear  deformation,  and  mini  sandwich  beam. 

15.  NUMBER  OF  PAGES 

42 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

UNCLASSIFIED 

18.  SECURITY  CLASSIFICATION 

OF  THIS  PAGE 

UNCLASSIFIED 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

UNCLASSIFIED 

20.  LIMITATION  OF  ABSTRACT 

Ul. 

NSN  7540-01-280-5500 


COMPUTER  GENERATED 


Standard  Form  298  {Rev.  2-89) 
Prescribed  by  ANSI  Srd  239-18 
2»t(H 


TABLE  OF  CONTENTS 


Introduction . . .  1 

1 .  Effect  of  Fiber  Cross-Section  on  Stiffness  and  Strength  of  Laminated  Composites... 2 

Micromechanics  Model  for  Elliptic  Fibers . 2 

Determination  of  the  Shear  Moduli  G12  and  G13 . 8 

Determination  of  the  Shear  Modulus  G23 . 10 

Stress  Concentration  Factor  at  the  Fiber/Matrix  Interface . 11 

Numerical  Results . . . 12 

Laminate  Analysis . 13 

Inplane  Elastic  Constants . 14 

Elastic  Constants  Under  Bending . 14 

Numerical  Results . 15 

2.  Failure  Criteria  for  Composites  Under  Combined  Load 

Experimental . 19 

Materials . 19 

Fabrication  of  Sandwich  Panels . 20 

Compression  Testing  of  Sandwich  Panels . 20 

Data  Reduction . 20 

Results  and  Discussion . 22 

3.  Stress  Analysis  Models  for  Thick  Composites .  27 

Modified  Shear  Deformation  Theory . 27 

Cylindrical  Bending  Under  Uniform  Lateral  Load . 31 

Numerical  Results . 35 

Discussion . 36 


INTRODUCTION 


This  report  covers  three  topics  relative  to  mechanics  issues  in  advanced  composites: 
effect  of  fiber  cross-section  on  stiffness  and  strength  of  laminated  composites;  failure 
criteria  for  composites  under  combined  loading;  and  stress  analysis  models  for  thick 
composites.  The  work  was  performed  for  WL/MLBM,  Wright-Patterson,  OH  45433-6533 
under  Contract  F33615-89-C-5609  with  Universal  Technology  Corporation.  Dr.  Richard 
Hall  was  the  project  monitor. 

In  the  area  of  fiber  cross-sectional  effect,  work  is  concentrated  on  looking  at  the  effect 
of  elliptical  fibers  on  the  stiffness  and  inplane  stress  concentrations  at  the  micromechanics 
level.  Calculated  ply  stiffness  properties  are  used  in  conjunction  with  classical  laminated 
plate  theory  to  determine  the  effect  of  fiber  geometry  on  laminate  inplane  and  bending 
stiffness.  The  micromechanics  model  is  based  on  a  strength-of-materials  approach  in 
conjunction  with  a  laminate  analogy.  All  nine  elastic  constants  of  the  unidirectional 
composite  are  determined.  Numerical  results  show  that  fiber  cross-sectional  geometry  has 
only  a  modest  influence  on  laminate  e’astic  constants. 

Failure  criteria  for  filament  dominated  laminates  subjected  to  compression  loading  is 
given  consideration  under  the  second  topic.  Of  particular  interest  is  the  interaction  between 
longitudinal  compression  loading  and  transverse  loading  (tension  and  compression).  Data 
obtained  in  conjunction  with  a  mini-sandwich  beam  are  used  to  generate  failure  stresses. 
Classical  lamination  theory,  including  thermal  residual  stresses  due  to  cure,  is  utilized  for 
determining  ply  stresses.  Residual  stresses  are  based  on  initial  elastic  properties,  while 
mechanical  stresses  are  based  on  secant  modulus  at  failure  (compression  stress-strain  curve 
for  both  unidirectional  and  laminated  composites  is  nonlinear  at  high  strain  levels). 
Experimental  results  reveal  little  interaction  between  longitudinal  compression  and 
transverse  stresses  as  long  as  the  transverse  load  is  below  uniaxial  transverse  failure.  A 
technical  paper  on  this  work  was  presented  at  the  American  Society  for  Composites 
Conference  held  at  Penn  State  University,  October  13-15, 1992.  The  paper  also  appeared 
in  the  conference  proceedings  published  by  Technomic  Publishing  Co.,  Lancaster,  PA. 

Under  the  topic  of  thick  composites,  the  effect  of  shear  deformation  and  transverse 
normal  stress  on  the  cylindrical  bending  of  laminated,  anisotropic  plates  subjected  to  a 
uniform  lateral  load  is  investigated.  Field  equations  are  based  on  Reissner’s  principle  in  a 
modified  form  which  requires  assumed  interlaminar  shear  and  normal  stress  distributions  in 
addition  to  the  inplane  displacements.  Closed  form  solutions  are  presented  for  laminates 
with  simply-supported  boundaries.  Numerical  results  indicate  that  the  effect  of  transverse 
normal  stress  on  maximum  deflection  is  small,  even  for  relatively  thick  plates. 
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1.  EFFECT  OF  FIBER  CROSS-SECTION  ON  STIFFNESS  AND 
STRENGTH  OF  LAMINATED  COMPOSITES 


The  effect  of  fibers  with  elliptical  cross-section  shape  on  composite  stiffness  and  local 
strength  is  investigated.  Analysis  is  accomplished  at  both  the  micromechanics  and 
macromechanics  (laminate)  level.  At  the  micromechanics  level  a  strength-of-materials 
model  based  on  a  laminate  analogy  is  utilized  to  determine  the  nine  elastic  constants  of  a 
unidirectional  composite  containing  continuous  fibers  with  elliptical  cross-section  shape. 
The  stress  concentration  factor  at  the  fiber-matrix  interface  under  transverse  normal  stress 
and  inplane  shear  stress  is  also  determined.  The  effect  of  fiber  geometry  on  these  stiffness 
properties  of  laminates  under  inplane  and  bending  loads  are  then  investigated  using  the 
results  of  the  micromechanics  model  in  conjunction  with  classical  laminated  plate  theory. 
Stiffness  results  are  compared  using  a  [0°/±  45°]s  laminate  geometry  in  conjunction  with 
plies  containing  various  combinations  of  circular  and  elliptic  fibers. 

Micromechanics  Model  for  Elliptic  Fibers 

Consider  a  continuous  unidirectional  composite  containing  fibers  with  an  elliptical 
cross-section.  A  repeating  element  of  the  composite  cross-section  is  assumed  to  be  of  the 
form  shown  in  Fig.  1.  Because  of  symmetry  it  is  only  necessary  to  consider  an  upper 


z 


Figure  1.  Analytical  micromechanics  model. 
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quarter  of  this  cross-section  in  the  analysis.  This  upper  quarter  is  divided  into  a  series  of 
plates  (or  plies)  of  thickness  dz.  Using  a  strength-of-materials  micromechanics  approach 
the  stiffnesses  of  each  of  these  plies  are  determined  as  a  function  of  fiber  properties,  matrix 
properties,  fiber  volume  fraction,  and  fiber  geometry.  The  dimension  A  can  be  determined 
as  a  function  of  a/b  and  the  volume  fraction  of  fiber,  Vf,  as  follows 

A  =  k  [-  (R  +  1 )  +  +  1)2  +  (JL  -  4)  1  (1) 

where  R  =  a/b. 


Oxfp 


Figure  2.  Ply  with  biaxial  normal  stress  loading 


3 


A  laminate  is  then  constructed  from  these  plies  and  the  resulting  stiffnesses  determined 
by  integrating  the  ply  properties.  The  effective  elastic  constants  are  determined  in  the  same 
manner  as  for  a  laminated  plate  [  1  ].  This  procedure  is  referred  to  as  the  laminate  analogy 
micromechanics  model  for  determining  elastic  constants  of  a  unidirectional  composite. 
This  is  an  extension  of  a  model  proposed  by  Eckvali  |2]  for  transverse  modulus  and 
inplane  shear  modulus  of  a  unidirectional  composite  containing  circular  fibers  in  a  square 
array  [2],  A  similar  model  was  used  by  Chamis  [3 J  for  determining  the  fiber-matrix 
interface  stress  concentration  factor  under  transverse  tensile  (or  compression)  load  and 
inplane  shear  load. 


Figure  3.  Through-the-thickness  ply  loading. 


We  now  consider  the  strip  shown  in  Fig.  2  for  the  ply  analysis.  The  subscripts  f,  m, 
and  p  refer  to  fiber,  matrix,  and  ply,  respectively.  Elastic  moduli  and  Poisson's  ratios  are 
determined  by  assuming  the  volume  average  inplane  ply  strains  e*  and  ey  to  be  constant 
through-the-thickness,  i.e  we  assume 

exp  =  Ex  =  constant,  £yp  =  Ey  =  constant  (2) 
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Axial  and  transverse  stresses  are  assumed  to  be  applied  and  will  vary  from  ply-to-ply 
because  of  the  change  in  local  ply  stiffness.  Because  of  continuity  requirements  at  the 
fiber-matrix  interface 

^yfp  =  ^ymp  =  ^yp  13) 

Through-the-thickness  ply  loading  is  shown  in  Fig.  3.  The  total  load  produced  by  a/p  is 
assumed  to  vanish,  i.e.  there  is  no  normal  force  transmitted  through-the-thickness.  Thus, 

^zfpVfp  +  OzmpO'Vfp)  =  0  (4) 

where  Vfp  is  the  volume  fraction  of  fiber  in  the  ply,  i.e. 

Vfp  =  — — —  (5) 

(a  +  A) 


and 


yp=aVr^f^ 

Assuming  a  transversely  isotropic  fiber  and  an  isotropic  matrix  the  stress-strain 
relations  for  the  constituents  are  of  the  form 

t^xfp  =  C 1 1  fEx  ■*"  Ci2f  (Cyfp  +  Ezp) 

Gyp  =  Ci2f€$  +  C22f€yfp  +  C23f€z.p 
®zfp  =  Ci3f£x  f-23f€yfp  +  C22f€zp 

(6) 

*7xmp  =  CnmEx  "*■  +  £tp) 

Gyp  =  Cl  2m  (Ex  +  £zp)  +  C]  lmEymp 
^zmp  =  Cj2m  (Ex  "E  Eymp)  +  CnmEzp 

where  Qj  are  stiffnesses.  It  should  be  noted  that,  although  the  volume  average  transverse 
strain  £ypis  identical  from  ply  to  ply,  eyp  is  not  the  strain  of  the  fiber  or  the  matrix.  In  fact 

Eyp  =  £y  =  Eyfp V fp  +  Eymp(  1  '  Vfp)  (7) 
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In  addition 


°xp  ~  dxfpVfp  +  ^xmpO  '  Vfp)  (8) 

Combining  eqs.  (4),  (7),  and  (8),  with  the  constitutive  relations,  eqs.  (6),  we  obtain 
the  following  relationships 

p  _  C-22fCl  ]  m _ 

22p  (Cnf<l- Vfp)l 

=  l-P1  VfP  +  ^  (1  -Vfp)l,  ^  =  [p23f  Vf  +  ^2®.  (i  -Vfp)l 

'-22p  v-22f  '-11m  C22p  <~22f  V-Ilm 


Ql Ip  =  QllfVfp  +  Qllm(l  -  Vfp),  Q22p  =  Q22fVfp  +  Ql  lm(l  -  Vfp) 


Ql2p  =  Ql2fVfp  +  Ql2m(l  -  Vfp) 


Cup  =  [Qiip  +  c22p  &&)\  C33p  =  iQ22p  +  c22p  (~£)2] 


Cl3p  =  [Q,2p  +  C22p(— ~)(§f^)J 


where  Q;j  are  reduced  stiffnesses  for  plane  stress.  The  reduced  stiffnesses  for  the  ply  can 
now  be  determined  utilizing  the  results  from  eqs.  (9).  In  particular. 


C13p2 


2v<r23Ih2l 


Qiip  =  [Cup - ~H,  Q22p  =  [C22p-(?.-L-)(C22pz)(-~li)  ] 

^33p  ^-33p  ^22P 

Ql2p  =  [C12p  -  (pr!-)(C22p2)(^)(^)l 

^33p  '-'22P  ^22P 


These  reduced  stiffnesses  are  now  averaged  through-the-thickness  to  produce  effective 
values  for  the  micromechanics  model,  i.e. 


Jrd^) 
r  b 

0 

D 


(i,j=  1,2) 
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This  integration  is  performed  numerical  using  100  strips  through-the-thickness.  Elastic- 
constants  can  now  be  determined  in  the  same  manner  as  in  laminated  plate  theory,  i.e. 


E^Qn-%4  E2  =  Q22-Qi22 


Q22 


Q 


Vi2  = 


11 


Ql2 

Q22 


(12) 


where  Ej,  E2,  and  v12  denote  the  modulus  parallel  to  the  fibers,  the  modulus  transverse  to 
the  fibers,  and  the  Poisson's  ratio  as  determined  from  a  uniaxial  tensile  test  parallel  to  the 
fibers  while  measuring  the  strain  transverse  to  the  fibers,  respectively. 

Recognizing  that  O/p  of  the  ply  vanishes  and  utilizing  the  ply  constitutive  relations,  we 
can  write 


Ct3p 

zp  =  ‘  c33p 


Ex 


^23p  0 

c33P  y 


(13) 


Integration  of  eq.  (13)  leads  to  an  effective  value  of  the  through-the-thickness  strain,  which 
can  be  written  in  the  form 


£7,  =  -  Hi36x  -  H23  £y 


(14) 


where 


Hi3  = 


_ 1 _ 

o  +£) 

b 


A) 

b  ^i3P  dz 
C33p  b 


(i=  1,2) 


We  now  consider  a  uniaxial  load  applied  parallel  to  the  fibers.  In  this  case  we  have 

By  =  -  v]2  6x,  e£  =  -vi3e2  (15) 

where,  obviously,  v13  denotes  the  Poisson's  ratio  as  determined  from  measuring  the  strain 
transverse  to  the  fibers  during  this  uniaxial  tensile  test  in  the  fiber  direction.  Combining 
eqs.  (14)  and  (15),  and  taking  into  account  eq.  (12),  we  find  that 

Vl3  =  H]3  -  H23  (16) 

Q22 
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In  a  similar  manner,  measuring  the  strain  through-the-thickness  during  a  uniaxial  tensile 
test  transverse  to  the  fibers  leads  to  the  Poisson’s  ratio  v23.  In  this  case  we  have 

e5  =  -vi2|^ey,  e?  =  -  v23  £y  (17) 

Combining  eqs.  (14)  and  (17),  and  taking  into  account  eq.  (12),  we  obtain  the  result 

v23  =  H23-H13S^  (18) 

Qn 

The  through-the-thickness  modulus,  E3,  can  be  calculated  by  rotating  the  model  in 
Fig.  1  by  90°,  i.e.  substituting  1/R  for  R  in  the  current  analysis  and  letting  E3  =  E2. 

Determination  of  the  Shear  Moduli  Gn  and  Gjq 

For  determination  of  the  inplane  shear  modulus  G12  we  consider  the  ply  loading 
shown  in  Fig.  4.  The  effective  shear  strain,  yxyp  is  assumed  to  be  the  same  within  each 
ply,  i.e. 


Fig  4.  Inplane  shear  loading  at  ply  level. 


8 


Yxyp  =  Txy  =  constant  (19) 

where  Hy  is  the  macroscopic  shear  strain  relative  to  the  x-y  plane. 

Continuity  of  the  shear  stress  at  the  fiber  matrix  interface  requires  tl> 


^xyf  —  'Exym  —  ^xyp 


(20) 


and  macroscopic  stress-strain  relation  is  of  the  form 


where 


(21) 


x 


xy  - 


_ 1 _ 

<1+t> 


Txyp 


dz_ 

b 


The  shear  strain  at  the  ply  level  can  be  shown  to  be  the  following  function  of  the  fiber  and 
matrix  shear  strains 


Yxyp  —  Txy  —  YxyfVfp  +  YxymO  '  Vfp) 


(22) 


The  constitutive  relations  for  the  fiber  and  matrix  are  of  the  form 


^xyp 


_  T’xyp 


Yxyf  ~  r  «  Yxym  -  p 
Gi2f  Gm 


(23) 


Combining  eqs.  (22)  and  (23),  we  obtain  the  ply  shear  modulus 


r  _  x*yp  _ 
Ci2p  =  — r—  = 


Gi2>Cir 


f  [Gi2f(l  -  Vfp)  +  GmVfp] 


(24) 


The  shear  modulus  G12  can  be  obtained  by  combining  eqs.  (21)  and  (24)  with  the  result 
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'12 


j-f 

+  A)  ;o 


d+^-) 
b 


Gi2iGr 


dz 


(l+£) 


[ G 1 2f(  1  -  Vfp)  +  GnjVfp]  b 


(25) 


The  shear  modulus  G13  can  be  obtained  in  the  same  manner  as  E3,  i.e.  Fig.  1  is 
rotated  90°  and  the  same  procedure  implemented  as  for  G12. 
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For  the  transverse  shear  modulus,  G23,  we  consider  the  through-the-ihickness  ply 
loading  as  shown  in  Fig.  5.  The  shear  strain  at  the  ply  level  is  assumed  to  be  constant,  i.e. 

Yyzp  =  Yyzf  =  Yyzm  (26) 

and  the  effective  ply  shear  stress  is  assumed  to  be  constant  through-out  the  thickness  of  the 
plate.  Thus 


^yzp  ~  %z 


(27) 


A  simple  force  balance  leads  to  the  result 

^yzp  =  ^yz  =  ^yzfVfp  +  tyzjn  ( 1  Vfp)  (28) 

This  relationship  does  not  maintain  continuity  of  Tyz  at  the  Fiber/matrix  interface  at  the  ply 


tyz 


Figure  5.  Through-the-thickness  shear  loading  at  ply  level 
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level.  However,  if  we  consider  a  very  steep  stress  gradient  near  the  fiber/matrix  interface, 
then  we  can  accept  eq.  (28)  on  an  average  stress  basis.  The  constitutive  relations  are  of  the 
form 


TyZ  —  G23pYyzp>  ^yzf  —  G23fYyzp*  ^yzm  —  GmYyzp  (29) 

Substituting  eq.  (29)  into  eq.  (28),  we  obtain  a  simple  rule-of-mixtures  relationship. 


G23p  =  G23fVfp  +  Gm(l-Vfp)  (30) 

The  effective  laminate  shear  strain  is  obtained  by  integrating  Yyzp  through-the-thickness 
with  the  result 


K  l+A. 

_  Tyz  _  Tyz  I  b  1  dz 

°23  (1  +  A)  Jo  G23>  b 

b 


(31) 


which  immediately  leads  to  the  relationship 


G23  = 


<1+f> 


r— l— 

J0  [G23fVfp  +  Gm  (1  -  Vfp)] 


dz 

b 


(32) 


We  have  now  determined  all  9  elastic  constants  for  the  repeating  elliptic  fiber  element. 

Stress  concentration  Factor  at  the  Fiber/Matrix  Interface 

The  stress  concentration  factor  (SCF)  at  the  fiber/matrix  interface  is  of  interest  for  both 
transverse  and  inplane  shear  loading.  For  these  loadings  the  maximum  stress  will  occur  at 
z  =  0,  y  =  b.  In  the  scheme  under  consideration  this  represents  the  first  element  in  the 
integration  procedure.  Thus,  for  transverse  loading 

<*yp(l)  =  Ql2p(l)e$  +  Q22p(l)Ey  (33) 
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where  (1)  denotes  the  first  element  in  the  integration  scheme.  Under  transverse  loading  we 
have 


(34) 


Utilizing  eqs.  (33)  and  (34)  in  conjunction  with  the  definition  of  SCF,  and  taking  into 
account  eq.  (17),  we  obtain  the  following  result 


SCF  = 


[Q22p(D-V,2|2-QI2P(l)] 
_ Ei _ 

E2 


(35) 


A  similar  procedure  in  conjunction  with  inplane  shear  loading  leads  to  the  following  SCF 
relationship 


SCF  = 


^xyp(  1)  _  Gi2p(l) 

TXy  G12 


(36) 


Numerical  Results 

In  the  results  presented  here  the  fiber  is  assumed  to  be  transversely  isotropic  and  the 
matrix  is  assumed  to  be  isotropic.  We  now  consider  the  following  constituent  properties 

Eif  =  35xl06  psi,  E2f  =  3xl06  psi,  Gi2f  =  5xl06  psi,  V]2f  =  0.25 
v23f  =  0.55,  Eif  =  35xl06  psi,  Em  =  5X105  psi,  vm  =  0.35 

Numerical  results  are  shown  in  Table  1  for  circular  fibers  and  for  elliptic  fibers  with 
various  a/b  ratios.  The  fiber  volume  fraction  is  assumed  to  be  0.6.  The  fiber  properties 
given  by  eq.  (37)  are  typical  of  state-of-the-art  graphite  fibers  and  the  matrix  properties  are 
those  of  high  performance  epoxy  resins. 

The  results  seem  in  Table  appear  to  be  reasonable.  There  is,  however,  an 
inconsistency  relative  to  G23.  In  particular  the  same  value  should  be  obtained  for  R  and 
1/R  aspect  ratios  (e.g.  R  =  2  and  0.5  should  yield  the  same  values).  This  is  not  the  case 
with  regard  to  the  results  shown  in  Table  1.  These  values  may,  however,  be  considered  as 
upper  and  lower  bounds. 

As  one  would  anticipate,  elliptical  fibers  with  R  >  1  increase  the  transverse 
modulus  and  inplane  shear  modulus,  while  the  same  moduli  are  reduced  for  R  <  1.  The 
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percentage  of  these  changes  are  not  great.  It  should  be  noted,  however,  that  the  low 
transverse  modulus  and  low  inplane  shear  modulus  of  the  fiber  contributes  to  the  smaller 
increase.  In  particular,  the  ratio  of  fiber  to  matrix  property  in  this  case  is  not  significantly 
large.  For  a  high  modulus  isotropic  fiber,  such  as  boron,  the  effect  of  fiber  cross-section 
shape  may  be  more  significant. 

TABLE  1  Micromechanics  Results  (Vf=  0.6) 


R=1 

R=2 

R=4 

R=l/2 

R=l/4 

Ex  (msi) 

21.190 

21.190 

21.190 

21.190 

21.190 

E2  (msi) 

1.452 

1.507 

1.562 

1.406 

1.375 

E3  (msi) 

1.452 

1.406 

1.375 

1.507 

1.562 

V12 

0.289 

0.286 

0.282 

0.293 

0.296 

V13 

0.289 

0.293 

0.296 

0.286 

0.282 

V23 

0.398 

0.413 

0.428 

0.383 

0.372 

G12  (msi) 

0.637 

0.728 

0.850 

0.577 

0.540 

G13  (msi) 

0.637 

0.577 

0.540 

0.728 

0.850 

O 

3 

'w' 

0.506 

0.480 

0.461 

0.535 

0.562 

SCF  (g2) 

1.456 

1.540 

1.621 

1.384 

1.332 

SCF  (X12) 

1.836 

2.085 

2.400 

1.663 

1.553 

It  should  also  be  noted  that  the  stress  concentration  factors  at  the  fiber  matrix 
interface  are  increased  for  R  >  1  and  decreased  for  R  <  1.  Thus,  it  is  anticipated  that  the 
transverse  strength  and  inplane  shear  strength  would  be  effected  in  a  corresponding 
fashion. 

Laminate  Analysis 

We  now  consider  the  effect  of  fiber  cross-sectional  shape  on  the  inplane  and  bending 
stiffnesses  of  laminates.  This  is  accomplished  by  combining  the  ply  elastic  properties 
produced  for  various  fiber  cross-sectional  shapes  (as  given  by  Table  1)  with  classical 
laminated  plate  theory. 
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Inplane  elastic  constants  for  symmetric  laminates  can  be  determined  from  classical 
laminated  plate  theory  by  considering  the  laminate  inplane  constitutive  relations  in  the  form 


£i  =  Qij£?  (i>  j  =  1.2,  6)  (38) 

where 


h/2 

rh/2 

C?°dz,  Qij  = 

h/2 

1  h/2 

and  the  superscript  k  denotes  the  kth  layer.  Then  the  inplane  effective  laminate  elastic 
constants  for  balanced  laminates  (same  number  of  angle-ply  layers  at  ±  0)  take  the  form 


Ex  =  Qii-2^-,  Ey  =  Q22-S12 


Q22 


Q11 


Gxy  —  Q66>  vxy  “  ^ 


Ql2 


xy 


Q22 


(39) 


Balanced  laminates  produce  orthotropic  inplane  elastic  constants  (no  shear  coupling).  This 
is  due  to  the  fact  that  the  effective  laminate  stiffnesses  are  independent  of  stacking 
sequence. 

Elastic  Constants  Under  Bending 


Under  bending  loads,  the  constitutive  relations  for  a  symmetric  laminate  are  of  the  form 


where 


s 

II 

0 

jj 

(bj  = 

1,2,  6) 

f  h/2 

/•h/2 

I  o^z  dz, 

DM  = 

J  Q-fz2  dz 

J- h/2 

* 

/  h/2 

(40) 


and  Kj  are  plate  curvatures.  In  order  to  obtain  effective  elastic  constants,  we  rewrite  eq. 


(40)  in  the  form 


Mi  -  J^QijKj  (4 1 ) 

where  the  superscript  b  denotes  bending  properties.  The  effective  bending  stiffnesses  are 
obviously  given  by 


Ob  =  11 D 

Q.J  h3  ■ 


Inverting  eq.  (41),  we  can  write  the  constitutive  relations  in  terms  of  compliances. 


Ki 


(42) 


where  S®  denote  effective  bending  coampliances  of  the  laminate. 

Since  the  bending  stiffnesses,  Djj,  are  a  function  of  stacking  sequence,  the  effective 
stiffnesses  for  balanced  laminates  will  not  be  orthotropic.  Thus,  shear  coupling  (in  the 
form  of  bending-twisting  coupling)  will  be  present  for  the  general  case  of  symmetric 
laminates  subjected  to  bending  loads.  Effective  bending  elastic  constants  can  be  determined 
from  the  compliances  in  eq.  (42)  in  the  usual  manner,  i.e. 


(43) 


where  i)*  and  are  bending  shear  coupling  coefficients. 


Numerical  Results 


Elastic  constants  under  inplane  and  bending  loads  are  shown  in  Table  2  for  a 
[0°2/±45°]s  laminates.  Ply  properties  are  taken  from  Table  1.  Four  different  combinations 
of  fiber  geometries  are  considered  as  follows: 


#1  All  plies  with  R  =  1 

#2  0°  plies  R  =  2,  45°  plies  R  =  1 
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#3  0°  plies  R  =  1/2,  45°  plies  R  =  1 

#4  All  plies  R  =  2 
#5  All  plies  R  =  1/2 


TABLE  2  Laminate  Elastic  Constants,  [0o2/±45o]s 


#i 

#2 

#3 

#4 

#5 

Ex  (Msi) 

11.900 

11.910 

11.890 

12.040 

11.800 

Ey  (Msi) 

3.419 

3.445 

3.398 

3.517 

3.349 

vxy 

0.711 

0.708 

0.714 

0.689 

0.727 

Gxy  (Msi) 

3.059 

3.104 

3.029 

3.109 

3.025 

Exb  (Msi) 

18.910 

18.910 

18.900 

18.940 

18.890 

Eyb  (Msi) 

1.883 

1.940 

1.836 

1.956 

1.825 

vxyb 

0.444 

0.441 

0.448 

0.434 

0.453 

Gxyb  (Msi) 

1.136 

1.218 

1.081 

1.220 

1.079 

^lxb 

-0.208 

-0.197 

-0.216 

-0.199 

-0.215 

V! 

a* 

-0.358 

-0.336 

-0.374 

-0.335 

-0.375 

As  easily  seen  from  Table  2,  fiber  geometry  has  little  effect  on  either  inplane  or 
bending  elastic  constants. 
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2.  FAILURE  CRITERIA  FOR  COMPOSITES 
UNDER  COMBINED  LOADING 


The  uniaxial  compression  strength  of  a  0°  composite  has  long  been  a  property  of 
significant  interest  to  both  materials  engineers  and  designers.  In  the  present  paper, 
however,  the  interaction  between  longitudinal  compression  and  transverse  loading  of  a  0° 
composite  is  investigated.  This  is  accomplished  by  subjecting  a  series  of  laminates  with 
different  stacking  geometries  to  inplane  compression  loading. 

Laminate  geometries  were  chosen  to  produce  a  number  of  longitudinal  compression  to 
transverse  stress  ratios.  In  particular,  for  multidirectional  laminates  the  0°  plies  will  be 
under  biaxial  loading  which  creates  an  interesting  challenge  relative  to  failure  criteria.  For 
example,  consider  the  three  loading  conditions  shown  in  Fig.  6.  One  would  anticipate  a 
higher  compression  strength  in  the  presence  of  transverse  compression  loading  (Fig.  6b) 
compared  to  pure  compression  (Fig.  6a).  However,  if  the  transverse  load  is  tension  (Fig. 
6c),  one  would  anticipate  a  lower  longitudinal  compression  strength.  We  envision  that 
transverse  compression  would  provide  additional  stability  to  the  longitudinal  compression 
loading,  while  transverse  tension  would  have  the  opposite  effect.  In  addition  transverse 
tension  could  possibly  induce  longitudinal  splitting  followed  by  fiber  microbuckling.  If  the 
dominate  failure  mode,  however,  is  fiber  failure,  then  little  interaction  between  longitudinal 
and  transverse  load  may  be  anticipated  provided  the  transverse  loads  are  below  the  uniaxial 
transverse  strength. 

The  anticipated  difference  in  longitudinal  compression  strength  between  the  loadings 
shown  in  Fig.  6  are  not  predictable  with  a  number  of  conventional  lailure  criteria. 
Maximum  stress  and  maximum  strain  do  not  recognize  coupling  between  longitudinal  and 
transverse  normal  stresses.  For  criteria  which  recognize  coupling  between  stress 
components,  such  as  the  quadratic  interaction  criterion  introduced  by  Tsai  and  Wu  [1], 
there  are  also  problems  which  we  will  now  consider.  This  criterion  is  of  the  form 


OL2  |  CTT2  Plot  ,  tLT2 

Slt  Slc  Sjt  Sjc  2Y  Slt  Slc  Syr  Stc  Tlt2 


+  (■ 


1 


1 


Slt  Slc 


■)  cl  +  ( 


j _ L 

Srr  Stc 


)  aT  =  1 


(44) 


where  oL,  Oj,  and  xLT,  are  the  normal  stresses  parallel  to  the  fibers,  transverse  to  the 
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fibers,  and  inplane  shear  stress,  respectively.  The  strength  parameters  SLT,  SLC,  S-pp.  Sxc, 
and  Tlt,  denote  tensile  strength  and  compression  strength  parallel  to  the  fibers,  tensile  and 
compression  strength  transverse  to  the  fibers,  and  the  inplane  shear  strength,  respectively. 


Fig.  6  Biaxial  Loading 


The  linear  terms  and  the  interaction  term  between  aL  and  aT  govern  the  loading  conditions 
shown  in  Fig.  6b  and  6c.  In  particular,  different  results  between  the  loadings  in  Fig.  6b 
and  lc  depend  on  a  distinct  difference  between  transverse  tension  and  compression 
strengths  and/or  differences  between  longitudinal  tension  and  compression  strengths.  A 
recent  paper  by  Theocaris  [2]  discusses  this  behavior  in  some  detail.  The  validity  of  eq. 
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(44)  in  conjunction  with  the  combined  loading  shown  in  Fig.  6  will  be  considered. 

Experimental 

Experimental  data  has  been  previously  obtained  by  the  University  of  Dayton  Research 
Institute  under  contract  F33  615-92-C-5618  with  the  Materials  Directorate  of  the  Air  Force 
Wright  Laboratory.  A  mini-sandwich  specimen  [3]  consisting  of  composite  face  sheets 
bonded  to  a  neat  resin  core,  as  illustrated  in  Fig.  7,  was  utilized  in  obtaining  this  data.  As 
opposed  to  the  original  sandwich  beam,  which  incorporated  a  honeycomb  core,  the  mini¬ 
sandwich  specimen  is  small,  relatively  inexpensive  to  fabricate,  and  uses  a  smaller  volume 
of  material. 


Fig.7  Mini-Sandwich  Beam  Specimen 


The  mini-sandwich  specimen  is  loaded  in  a  conventional  IITRI  test  Fixture  (ASTM  D- 
3410)  and  the  data  reduced  using  classical  laminated  plate  theory.  The  materials, 
fabrication  procedures,  and  test  procedures  utilized  by  the  University  of  Dayton  Research 
Institute  are  described  in  the  following  paragraphs. 

Materials 

The  cores  of  the  miniature  sandwich  specimens  were  fabricated  from  a  neat  epoxy  resin 
(Epon  828  from  Shell  Chemical  Co.),  cured  with  Jeffamine  D-230  (a  polyetheramine  from 
Texaco  Inc.).  The  composite  skins  were  fabricated  from  a  graphite/epoxy  AS4/350'. -6 
prepreg  (Hercules,  Inc.). 
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Fabrication  of  Sandwich  Panels 


The  miniature  sandwich  panels  were  fabricated  in  a  two-step  process.  In  the  first  step, 
the  following  symmetric  laminate  panels,  each  12  in.  x  12  in.,  were  first  fabricated  :  |()]2, 
[02/±30]s,  [±45/0/90] s,  [0/45/90/-45]s,  and  [0/±60]s.  Two  plates,  each  10  in.  x  6  in. 
were  cut  from  each  panel.  These  were  lightly  sanded  to  remove  any  residual  release  agent. 
A  mold  was  constructed  by  clamping  these  plates  across  a  uniform  Teflon  spacer  (3.18  mm 
thick)  along  three  sides  to  produce  a  rectangular  cavity  open  at  one  end.  Approximately 
125  g  of  Epon  828  were  mixed  with  32  phr  Jeffamine  D-230,  and  debulked  under  vacuum 
to  remove  entrapped  air.  The  liquid  resin  was  then  cast  in  this  mold  and  allowed  to  cure  at 
room  temperature  for  7  days.  The  product  was  a  symmetric  miniature  sandwich  panel  with 
excellent  bonding  between  composite  skin  and  resin  core.  In  order  to  measure  transverse 
compression  strength,  a  [90°]24  solid  panel  (no  core)  was  also  fabricated. 


j wL 


ompression  Testing  of 


Sandwich.  Panels 


After  assuring  laminate  quality  by  ultrasonic  C-scans,  compression  tests  were 
performed  on  specimens  cut  from  these  sandwich  panels.  Specimen  sizes  were  in 
accordance  with  ASTM  standards  for  conventional  laminates,  and  tests  were  also 
conducted  under  similar  conditions.  Axial  compression  testing  was  performed  in  an  IITRI 
fixture  on  specimens  127  mm  long,  and  6.4  mm  wide.  Specimens  were  gripped  through 
bonded  tapered  glass/epoxy  tabs  with  a  gage  length  of  12.7  mm.  Strain  gages  were 
bonded  to  both  faces  of  each  specimen  to  monitor  any  bending  or  buckling  resulting  from 
eccentric  loading.  An  additional  strain  gage  was  bonded  to  the  side  of  selected  specimens 
to  simultaneously  measure  axial  strain  in  the  core.  Tensile  tests  were  also  performed  on 
specimens  from  these  panels.  Failed  specimens  were  examined  in  an  optical  microscope 
and  also  in  a  JEOL  JSM  840  scanning  electron  microscope. 

Datfl-B.gdutiiQn 


Data  was  reduced  using  classical  laminated  plate  theory.  This  procedure  involved  initial 
use  of  macroscopic  constitutive  relations,  which  are  of  the  form 


'  N, 

An 

Ai2 

£i° 

NiT  " 

.  n2 

.  Ai2 

A22  . 

-  £2°  - 

.  N2T  _ 
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where  as  in  eq.  (44)  the  subscripts  L  and  T  refer  to  a  coordinate  system  parallel  and 
transverse  to  the  fibers.  Thermal  expansion  coefficients  in  the  longitudinal  and  transverse 
directions  are  denoted  by  a l  and  <Xt,  respectively.  The  value  of  AT  is  chosen  as  the 
difference  between  room  temperature  and  cure  temperature.  For  mechanical  stresses  the  0° 
secant  modulus,  Ej,  at  failure  is  utilized,  while  the  other  elastic  properties  are  assumed  to 
be  constant.  Since  the  core  is  not  co-cured  with  the  face  sheets,  thermal  stresses  are  based 
on  cure  of  the  face  sheets  independent  of  the  core. 

For  mechanical  stresses  the  longitudinal  compression  secant  modulus  of  the 
unidirectional  ply  and  the  compression  secant  modulus  of  the  core  at  failure  are  used  due  to 
the  nonlinear  response  of  these  properties.  Other  elastic  properties  are  assumed  to  be 
constant.  The  Poisson's  ratio  of  the  core  is  assumed  to  be  v  =  0.35.  In  the  determination 
of  mechanical  stresses,  NjT  and  N2  (i.e.  the  load  is  applied  in  the  x\  direction)  are  allowed 
to  vanish.  Typical  stress-strain  curves  in  compression  for  the  core  and  a  unidirectional 
composite  are  shown  in  Figs  8  and  9,  respectively.  Because  most  of  the  load  is  carried  by 
the  unidirectional  face  sheets,  the  nonlinear  response  seen  in  Fig.  9  is  due  primarily  to  the 


compression  response  of  the  0°  composite. 


Fig.  8  Compressive  Stress-Strain  Curve  for  Neat  Epon  828  Resin 


Results  and  Discussion 

Experimental  results  are  shown  in  Table  3.  All  data  are  extrapolated  to  60%  fiber 
volume.  The  compression  strength  of  the  face  sheet  in  the  x  and  y  directions  are  denoted 
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by  Sxf,  and  Syf,  respectively,  while  the  0°  ply  stresses  at  failure  (including  thermal  residual 


Fig.  9  Stress-Strain  Curve  for  0°  Mini-Sandwich  Beam 


stresses)  in  the  longitudinal  and  transverse  directions  are  denoted  by  <7l  and  Oj, 
respectively.  Longitudinal  and  transverse  residual  stresses  of  the  0°  plies  are  denoted  by 
Olt  and  <Jtt,  respectively.  Note  that  due  to  a  mismatch  between  the  effective  Poisson's 
ratio,  vxy,  of  the  face  sheet  and  the  core,  the  face  sheets  will  be  under  biaxial  inplane 
normal  stresses  at  the  macroscopic  (laminate)  level.  In  addition,  there  are  no  inplane  shear 
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stresses  present  in  the  0°  plies.  Thus,  the  unidirectional  plies  are  only  subjected  to  biaxial 
normal  stresses  inplane.  Results  for  the  [0/45/90/-45]s  face  sheets  are  not  shown,  as  the 
initial  failure  was  due  to  ffee-edge  delamination  which  reduced  the  apparent  compression 
strength  considerably. 


TABLE  3  Experimental  Data,  psi 


Laminate 

Sxf 

Syf 

KM 

Gj 

°lt 

axT 

[90o]24* 

-32,800 

0 

0 

-32,800 

0 

0 

[0°]2 

-281,550 

765 

-281,550 

765 

0 

0 

[0°/90°]s 

-141,750 

7,964 

-262,680 

-213 

-4,780 

4,780 

[±45°/0790°]s 

-103,745 

454 

-260,680 

5,027 

-4,780 

-4,780 

[0y±30°ls 

-135,145 

-3,241 

-166,340 

8,255 

4,564 

1,966 

[0°£60°]s 

-97,208 

740 

-290,612 

5,215 

-4,780 

4,780 

♦Solid  laminate  tested  for  transverse  compression  strength 


Experimental  results  from  Table  3  were  compared  to  the  Tsai-Wu  failure  criterion,  eq. 
(44).  For  the  ply  loading  under  consideration,  eq.  (44)  reduces  to  the  form 


ctL2  l  Or2 _ <?l  C7j _ 

Slt  Slc  Stt  Stc  2V  Slt  Slc  ^tt  Stc 


1 

Slc 


■)  ol  +  ( 


1 

Sn 


1 

Stc 


■)  ot  =  1 


(48) 


The  following  ply  strength  parameters  were  used  in  comparing  theory  and  experiment: 

Slt  =  250  ksi,  Slc  =  282  ksi,  Stt  =  7.5  ksi,  Stc  =  33  ksi  (49) 


A  comparison  between  eq.  (48)  and  the  experimental  data  is  shown  in  Fig.  10.  Except 
for  the  results  of  the  [0°2/±30°ls  laminate,  the  longitudinal  compression  strength  appears  to 
be  relatively  insensitive  to  transverse  load.  In  the  case  of  the  [O02/±3()°]s  laminate  the 
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transverse  tensile  stress,  ax,  exceeds  the  transverse  tensile  strength,  Srr.  When  transverse 
tensile  failure  occurs  first,  the  fibers  lose  support  from  the  matrix  and  fiber  buckling 
becomes  the  failure  mode.  For  the  other  laminates  the  predominant  failure  mode  appears  to 
be  shear  failure  of  the  fiber  as  previously  described  by  Crasto  and  Kim  [2].  Thus  if  Cj  < 
Stt,  then  little  or  no  interaction  exists  between  longitudinal  compression  strength  and 
transverse  tensile  loading.  One  caution  must  be  raised  here,  however,  because  the  laminate 
geometries  presented  here  do  not  induce  transverse  compression  in  the  0°  plies  (when 
thermal  residual  stresses  are  included),  other  than  the  very  small  transverse  compression 
stress  present  in  the  [0°/90°]s  laminate. 


Fig.  10  Comparison  Between  Eq.  (48)  and  Data  in  Table  3 
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Thus  a  meaningful  conclusion  cannot  be  made  concerning  the  interaction  between 
longitudinal  compression  and  transverse  compression.  Such  data  is  very  difficult  to  obtain 
in  flat  laminates,  leaving  combined  axial  compression  and  external  pressure  loading  of 
tubular  specimens  as  the  only  likely  means  of  generating  this  information. 

It  is  also  important  to  note  that  0°  compression  strength  is  very  high  and  carries  over  to 
laminate  in-situ  unidirectional  plies. 
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3.  STRESS  ANALYSIS  MODELS  FOR  THICK  COMPOSITES 


It  is  well  recognized  that  shear  deformation  can  be  more  significant  in  laminated, 
anisotropic  plates  than  in  isotropic,  homogeneous  plates.  This  is  due  to  the  large  ratios  of 
effective  inplane  tensile  moduli  to  through-the -thickness  shear  moduli.  Shear  deformation 
theories  based  on  Mindlin's  assumptions  [4]  have  been  developed  for  laminated, 
anisotropic  plates  [5,  6].  A  theory  has  also  been  developed  [7 j  for  the  bending  of 
laminated,  anisotropic  plates  which  includes  the  effects  of  transverse  shear  deformation  in  a 
manner  similar  to  Reissner's  theory  for  homogeneous,  isotropic  plates  [8].  This  theory  is 
based  on  Reissner's  principle  [9]  in  a  modified  form  which  requires  assumed  distributions 
of  inteilaminar  shear  stresses  in  addition  to  the  inplane  displacements.  This  theory  has 
been  modified  [10]  to  include  the  effect  of  transverse  normal  stress  without  increasing  the 
number  of  bending  equations  above  those  generated  in  classical  shear  deformation  theory. 
In  this  theory  the  interlaminar  normal  stress  distribution  is  assumed,  in  addition  to  the 
interlaminar  shear  stress  distributions.  Numerical  results  for  cylindrical  bending  of 
laminated,  anisotropic  plates  subjected  to  a  half  sine  wave  lateral  load  showed  that 
improved  results  for  bending  deflections  can  be  obtained  by  including  transverse  normal 
stress  effects  as  well  as  transverse  shear  deformation  effects.  These  numerical  results  were 
obtained  for  simply-supported  boundary  conditions  and  results  compared  to  exact  solutions 
from  classical  theory  of  elasticity. 

In  the  present  effort  the  effect  of  transverse  shear  and  normal  stresses  on  the  bending 
of  symmetrically  laminated,  anisotropic  plates  subjected  to  cylindrical  bending  under 
uniform  lateral  load  is  investigated.  Closed  form  solutions  are  developed  for  simply- 
supported  boundary  conditions. 

Modified  Shear  Deformation  Theory 

For  continuity  purposes,  some  of  the  details  in  the  development  of  the  field  equations 
developed  by  Whitney  and  Rose  [10]  will  be  outlined  here.  Consider  a  rectangular 
laminated  plate  of  thickness  h  with  a  symmetric  stacking  sequence  relative  to  the  plate 
midplane.  Each  ply  is  constructed  of  an  orthotropic  material  with  the  principle  axis  of 
orthotropy  oriented  at  an  arbitrary  angle,  0,  relative  to  the  x-axis  of  the  plate.  The  inplane 
displacements  are  assumed  to  be  of  the  form 

u  =  u°(x,  y)  +  z  \|/X(x,  y),  v  =  v°(x,  y)  +  z  yy(x,  y)  (50) 
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while  the  plate  deflection,  w,  is  expressed  in  terms  of  the  weighted  displacement 


w* 


w(x,y,z)[l  -4(£)2]dz 


(51) 


The  ply  constitutive  equations  are  of  the  form 

Oj  =  CikEk,  £j  —  SjkOfc  (52) 

where  C&  and  Sjk  are  the  stiffnesses  and  compliances,  respectively.  We  now  eliminate  the 
interlaminar  normal  strain  £3  by  utilizing  the  first  of  the  constitutive  relations,  eq.  (52),  in 
the  form 


Oi  =  Qik£k  +  ^  ct3  (i,  k  =  1,  2,  6)  (53) 

C33 

where  Qij  are,  as  previously,  the  reduced  stiffnesses  for  plane  stress  given  by  the 
relationship 


Qik  =  Cik-^#2- 

<-33 

Reissner's  principle  is  written  in  terms  of  the  variational  equation  61 1  =  0.  Omitting 
boundary  and  surface  traction  terms  for  simplicity,  the  variational  function  811  can  be 
written  in  the  modified  form  [10] 

8n  =  J  { [<Ji  -  (QijEj  +  ^  03)]  5e;  +  (£j  -  Skm0m)  SOj 

+  (ai,i  +  06,2  +  05,5)  8u  +  (06,1  +  02,2  +  04,4)  8v 


i,  j  =  1,  2,  6' 
k,  m  =  4,  5 


(54) 


(05,1  +  04,2  +  03,3)  5w  }  dV 

where  a  repeated  index  denotes  summation  and,  in  the  usual  manner,  0;  and  £j  denote  the 
stress  and  strain  components,  respectively.  Using  conventional  nomenclature,  the 
subscripts  1 ,  2,  and  3  here  denote,  x,  y,  and  z,  respectively  while  4,  5,  and  6  correspond 
to  yz,  xz,  and  xy,  respectively. 
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We  now  denote  normal  stresses,  shear  stresses,  normal  strains,  and  shear  strains  by 
<Ji,  Tij,  £j  and  Yij,  respectively.  In  order  to  implement  eq.  (54),  we  need  to  assume  a 
distribution  through  the  thickness  for  xxz,  tyz,  and  0Z.  As  in  the  original  work  by  Whitney 
and  Rose  [10],  we  assume  these  stress  components  to  be  of  the  same  form  as  utilized  by 
Reissner  [8]  for  homogeneous,  isotropic  plates,  namely 

tyz  =  [l-4(f)2)^  <55) 

Gz  =  ^  {  [3  -  4  (J)2]  (J)  (pi  -  p2)  +  (pi  +  p2)  } 

where  Q*  and  Qy  are  the  shear  force  resultants  defined  in  the  usual  manner 

fh/2 

(Qx j  Qy)  =  I  (^xz>  tyz)  dz 
J  h/2 


and 


Pi  =  0z(h/2),  p2  =  0z(-h/2) 


These  relationships  are  consistent  with  linear  inplane  stresses  through-the-thickness 
for  the  case  of  homogeneous  plates.  It  should  be  noted  that  the  incorporation  of  az  in  the 
constitutive  relations,  eq.  (53),  will  induce  some  nonlinearity  through-the-thickness  with 
regard  to  the  inplane  stresses.  This  may  be  of  little  consequence,  however,  in  obtaining 
macroscopic  response  of  the  laminate.  Detailed  local  interlaminar  stress  distributions 
require  integration  of  the  equations  of  equilibrium  in  conjunction  with  the  inplane  stresses 
as  determined  from  the  laminated  plate  theory. 

Substituting  eqs.  (50),  (51),  and  (55)  into  eq.  (54),  we  obtain  the  following  inplane 
constitutive  relations  in  terms  of  force  resultants  and  moments 


where 


Nj  —  Ay  £j  +  Lj  (pi  +  p2) 
M;  =  Dy  Kj  +  J;  (pi  -  p2) 


(i.  j  =  1,2,6) 


(56) 
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(Ni,  Mi)  = 


/•h/2 

r  h/2 

1  ai(l,z)dz, 

(Aij,  Djj)  =  Qij  (1,  z2)  dz 

J  h/2 

J-  h/2 

h/2 

rh/2 

Si3  dz,  Ji  =  i 

Si3  (u  )[3  -  4  (l)1]  z  dz 

C33  2 

h/2 

C33  n  1  V  J 

f-  h/2 

We  note  in  eq.  (56)  that  the  inplane  force  resultants  are  coupled  to  the  lateral  surface 
tractions.  Thus,  the  midplane  strains  do  not  vanish  in  the  case  of  the  bending  due  to  lateral 
load  when  the  effects  of  transverse  normal  stress  are  considered.  Reissner  [8]  did  not 
consider  midplane  strains  in  the  development  of  his  theory  for  homogeneous,  isotropic 
plates.  In  a  similar  manner  we  obtain  the  constitutive  relations  for  the  transverse  shear 
force  resultants,  Qx  and  Qy,  from  eq.  (54)  with  the  result 


[Qx" 

F55 

F45 ' 

Vx  +  W*,x' 

LQy. 

L  F45 

F44  . 

\j/y  +  w*,y 

where 


(57) 


f44  = - ^ - ,  f55  = - ^ - ,  f45  =  -  ■  - 

(R44  r55  -  R452)  (R44  R55  -  R452)  (R44  R55  -  R452) 


and 


[1-4(J)2]2  dz 


(i,  j  =4,  5) 


Here  Sjj  denotes  the  anisotropic  transverse  shear  compliances. 

The  equilibrium  equations  in  terms  of  force  and  moment  resultants  can  be  obtained 
from  the  last  three  terms  in  eq.  (5)  with  the  result 

NX,X  +  Nxy,y  =  0 

NXy>X  Ny,y  =  0 

Mx,x  +  Mxy,y  -  Qx  =  0  (58) 

MXy,x  +  My,y  *  Qy  =  0 
Qx.x  +  Qy,y  +  (Pi  -  P2)  =  0 


These  five  equations  could  be  expressed  in  terms  of  the  displacement  variables  u°,  vO, 
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yx,  Vy,  and  w*  by  utilizing  the  constitutive  relations,  eqs.  (56)  and  (57).  The  boundary 
conditions  are  those  of  classical  shear  deformation  theory  as  applied  to  laminated  plates  [6|. 

Cylindrical  Bending  Under  Uniform  Lateral  Load 

For  cylindrical  bending  we  consider  a  symmetric  laminate  where  the  length  in  the  y- 
direction  is  much  longer  than  the  x-direction,  as  illustrated  in  Fig.  11.  If  the  laminate  is 
subjected  to  a  lateral  load  which  is  independent  of  y,  then  the  resulting  deformations  will  be 
independent  of  y.  Such  cases  are  referred  to  as  cylindrical  bending. 


Fig.  11  Cylindrical  Bending 

The  displacements  in  eq.  (1)  and  (2)  now  take  the  form 

u°  =  U(x),  v«  *  V(x),  yx  =  4/x(x),  \]/y  =  vFy(x),  w*  =  W(x)  (59) 
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Denoting  differentiation  with  respect  to  x  by  a  prime,  the  first  two  equilibrium 
equations  in  (58)  reduce  to  the  following 

Nx'  =  0,  Nxy'  =  0  (60) 

w’hich  leads  to  constant  values  of  Nx  and  Nxy.  Combining  eq.  (60)  with  the  constitutive 
relations,  eqs.  (56)  and  (57),  and  substituting  the  results  into  the  last  three  equations  in 
(58),  we  obtain  the  following  bending  equations 

Du  Hy  -  f55  +  d16  ny  -  f45  -  f55  w  +  jx  (Pi'  -  P2')  =  o 

D,6  'Fx"  -  F45  'Fx  +  D66  VFy"  -  F44  *Fy  -  F45  W'  +  J6  (Pr  -  p2')  =  0  (61) 

F55  'Fx*  +  F45  'Fy*  +  F55  W"  +  (Pi  -  p2)  =  0 


For  a  uniform  load 


pi  =  po  =  constant,  p2  =  0  (62) 

Because  the  lateral  load  is  constant,  the  effect  of  transverse  normal  stress  vanishes  from  the 
first  two  equations  in  (61),  with  the  result 

Du  'Fx"  -  F55  Wx  +  D16  'Fy"  -  F45  'Fy  -  F55  W' 

Die  'Fx"  -  F45  'Fx  +  D66  'Fy'*  -  F44  'Fy  -  F45  W* 

F55  'Fx’  +  F45  'Fy'  +  F55  W"  +  po  =  0 

The  solution  to  eqs.  (63)  is  given  by 

\ifx  =  ci  +  c2x  +  x2  -  (A  cosh  +  B  sinh  x3 

2  Du  h  h  6Dn 


=  0 

=  0  (63) 


\i/y  =  Hi  (C3  -  +  A  cosh  +  B  sinh 

1  Du  h  n 


(64) 


P«_x2 


w  =  c4  -  Cix  -  X2  +  (6H3  -  x2)  x  -  (12H2  -  X2) 

2  6  24L)it 

+  ^2.  (A  sinh  X£  +  B  cosh  X~) 

X  h  h 
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where 


x  =  /(F44F55-F452)Diih2  Hj  _(F55Di6-F45Dii) 

V  F55(D„D66  -  D162)  ’  1  (F44F55  -  F452) 

H2  -(P44D11  -  F45D}6)  pj  _  (F55D16  -  F4sDn) 

2  (F44F55  -  F452)  ’  3  p55Dn 

The  width  of  the  plate  is  denoted  by  a  and  the  axis  system  is  placed  at  the  middle  such 
that  the  boundaries  are  located  at  x  =  ±  a/2  (see  Fig.  11).  At  the  laminate  level,  transverse 
normal  stress  will  only  effect  cases  in  which  force  and/or  moment  boundary  conditions  are 
required.  In  particular,  because  transverse  normal  stress  effects  do  not  appear  in  eqs.  (63), 
they  enter  the  problem  through  force  and/or  moment  boundary  conditions,  as  observed  in 
the  constitutive  relations,  eqs.  (56).  Thus,  for  example,  in  the  case  of  clamped  boundaries, 
which  involves  displacement  conditions  only,  transverse  normal  stress  effects  will  not  be 
present  at  the  laminate  level.  Transverse  normal  stress  will,  however,  have  an  effect  on  the 
ply  stresses  by  virtue  of  eq.  (53). 

We  now  consider  simply-supported  boundary  conditions.  For  this  case  w  will  be 
symmetric  about  x  =  0,  while  both  \|/x  and  \j/y  are  anti-symmetric.  This  leads  to  the  result 


Ci  =  C3  =  A  =  0  (65) 

For  simply-supported  boundaries,  we  require  that  at  x  =  ±  a/2 

w  =  Mx  =  Mxy  =  Nx  =  Nxy  =  0  (66) 

where,  from  the  constitutive  relations,  eqs.  (66) 

Mx  =  Dn\|/XiX  +  Di6\j/y>x  +  Jipo 

MXy  =  Dl6\jfXtX  +  D66\|tyjX  +  J6P0 


Combining  eqs.  (64)  and  (67),  we  find  that 
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c^_mlDl6  +  Du^-^lDn)  ^  B_  h  [H  Du  (JftDii  -  JjDiftj 

8Dn2  DnXcoshX^-  (DuD66  -  D]62) 

2h 

(68  "i 

c4  =  {[48(H2D„  +  HiD16  -  JiDu)  +  D„a2]Xa 

-  384h2H3Dn  [Hj  -  D11.-(J^PlI:. ilPli  ]  j - - 

(DnD66  -  Di62)  384X2Dn2 

The  maximum  deflection  will  occur  at  the  center  of  the  plate,  x  =  0,  and  can  be  written  in 
the  form 


Wmax  =  WCLPT(l  +  y) 


(69) 


where 


H  =  %H2  +  eiiH,-  J,)- 

a2  On 


384h2H3 

XV 


(Di iD66  -  Di62)  cosh 


and  wclpt  is  the  maximum  deflection  as  determined  from  classical  laminated  plate  theory 
with  shear  deformation  and  transverse  normal  stress  effects  neglected  and  is  of  the  form 

WCLFr  =  384D7T  P"  (70) 

Thus,  the  term  H  contains  the  effects  of  both  transverse  shear  deformation  and 
transverse  normal  stress.  We  note  that  the  classical  solution  is  independent  of  the  bending¬ 
twisting  coupling  stiffness  D^.  For  the  case  of  an  orthotropic  material  F45  =  Di6  =  0  and 
eq. (69) reduces  to 


wmax  =  wclpt  [1 


+  48(Pn  -  J1F55) 
5F55a2 


(71) 


We  should  also  note  that  midplane  strains  will  also  occur.  Combining  the  boundary 
conditions,  eq.  (66),  with  eq.  (60),  we  find  that 

Nx  =  Nxy  =  0  (72) 
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throughout  the  plate.  Substituting  eq.  (72)  into  the  first  of  the  constitutive  relations,  eq. 
(56),  we  obtain  the  midplane  strains 


„  o  _  (A6tLi  -  AjgLg) 

0 

(AuA66  '  Ai6^) 
o  _  (AnL6  -  Ai6L]) 

Yxy  — 

(AhA66  -  Aj62) 


(73) 


Numerical  Results 


We  now  consider  laminates  with  the  stacking  geometry  [0°/90°]s  and  [±45°JS  having  the 
following  unidirectional  properties 


El/Et  —  14,  Glt/Et  =  Gu/Ej  =  0.615,  Gt3/Ej  =  0.323, 
VLT  =  Vl3  =  0.3,  Vt3  =  0.55 


(74) 


where  the  subscript  L,  T,  and  3  denote  axes  parallel  to  the  fiber,  perpendicular  to  the  fiber, 
and  through-the-thickness  of  a  unidirectional  composite,  respectively.  In  addition  E;,  Gjj, 
and  Vjj  denote  the  Young's  modulus  in  the  ith  direction,  the  shear  modulus  relative  to  the  i-j 

TABLE  4  Maximum  Deflection 
Wmax/WCLPT 


[0°/90°]s 

[±45°; 

Is 

a/h 

WSN 

WSD 

WSN 

WSD 

5 

1.901 

1.923 

1.421 

1.442 

10 

1.225 

1.231 

1.107 

1.112 

20 

1.056 

1.058 

1.027 

1.028 

30 

1.025 

1.026 

1.012 

1.013 
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plane,  and  the  Poisson's  ratio  measured  by  contraction  in  the  jth  direction  during  a  uniaxial 
tensile  test  in  the  ith  direction,  respectively.  These  properties  represent  state-of-the-art 
unidirectional  graphite/epoxy  composites. 

Numerical  results  are  shown  in  Table  4  for  the  maximum  deflection.  Equation  (7 1 )  is 
applicable  to  the  [0°/90°]s  stacking  geometry,  while  eq.  (20)  is  required  for  the  anisotropic 
[±45°]s  laminate.  The  normalized  deflection,  wSN,  represents  eqs.  (6(0  and  (71)  with  both 
shear  deformation  and  normal  stress  effects  included,  while  wSD  represents  these  same 
equations  with  normal  stress  effects  neglected  (Jj  =  J6  =  0). 

Discussion 

A  cursory  examination  of  table  1  reveals  that  the  inclusion  of  transverse  normal  stress 
reduces  the  maximum  deflection  compared  to  shear  deformation.  This  is  consistent  with 
the  results  presented  by  Whitney  and  Rose  [10]  where  comparison  with  exact  elasticity 
solutions  showed  that  shear  deformation  alone  overestimated  the  maximum  deflection.  The 
inclusion  of  both  shear  deformation  and  transverse  normal  stress,  in  the  same  manner 
presented  here,  gave  almost  identical  results  to  exact  elasticity. 

The  effect  of  shear  deformation  and  transverse  normal  stress  is  more  severe  in  the 
[0°/90°]s  laminate  than  in  the  [±45°]s  stacking  geometry.  This  is  anticipated  as  the  [0°/90o]s 
laminate  has  a  higher  effective  inplane  modulus  to  interlaminar  shear  modulus  ratio  than  the 
[±45°]s  angle-ply  composite. 
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